Methods for nanosecond-scale time synchronization over a network

ABSTRACT

A method includes, at a first node: transmitting a first synchronization signal at a first time according to a first clock of the first node; back-coupling the first synchronization signal to generate a first self-receive signal; calculating a time-of-arrival of the first self-receive signal according to the first clock; and calculating a time-of-arrival of the second synchronization signal according to the first clock. The method also includes, at the second node: transmitting the second synchronization signal at a second time according to a second clock of the second node; back-coupling the second synchronization signal to generate a second self-receive signal; calculating a time-of-arrival of the second self-receive signal according to the second clock; and calculating a time-of-arrival of the first synchronization signal according to the second clock. The method S 100  further includes calculating a time bias and a propagation delay between the pair of nodes based on the time-of-arrivals.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims the benefit of U.S. Provisional Application No. 62/668,219, filed on 7 May 2018, and U.S. Patent Application No. 62/734,978, filed on 21 Sep. 2018, both of which are incorporated in their entireties by this reference.

TECHNICAL FIELD

This invention relates generally to the field of networking and digital communication and more specifically to a new and useful method for network time synchronization in the field of networking and digital communication.

BRIEF DESCRIPTION OF THE FIGURES

FIG. 1A is a flowchart representation of a method;

FIG. 1B is a flowchart representation of a first variation of the method;

FIG. 1C is a flowchart representation of a second variation of the method;

FIG. 2 is a schematic representation of a mesh network;

FIG. 3A is a schematic representation of a pair of nodes in the mesh network;

FIG. 3B is a schematic representation of a node in the mesh network;

FIG. 3C is a schematic representation of self-receive hardware of a node;

FIGS. 4A and 4B are schematic representations of synchronization slot structures;

FIG. 5A is a flowchart representation of transmission and reception of a synchronization signal; and

FIG. 5B is a flowchart representation of transmission and reception of a self-receive signal.

DESCRIPTION OF THE EMBODIMENTS

The following description of embodiments of the invention is not intended to limit the invention to these embodiments but rather to enable a person skilled in the art to make and use this invention. Variations, configurations, implementations, example implementations, and examples described herein are optional and are not exclusive to the variations, configurations, implementations, example implementations, and examples they describe. The invention described herein can include any and all permutations of these variations, configurations, implementations, example implementations, and examples.

1. METHODS

As shown in FIG. 1A, a method S100 for characterizing time bias and propagation delay between a pair of nodes includes, at a first node in the pair of nodes and during a first synchronization slot: transmitting a first synchronization signal at a first time according to a first clock of the first node in Block S110A; back-coupling the first synchronization signal to generate a first self-receive signal in Block S120A; receiving the first self-receive signal in Block S130A; calculating a time-of-arrival of the first self-receive signal according to the first clock in Block S140A; receiving a second synchronization signal from the second node in Block S150A; and calculating a time-of-arrival of the second synchronization signal according to the first clock in Block S160A. The method S100 also includes, at the second node and during the first synchronization slot: transmitting the second synchronization signal at a second time according to a second clock of the second node in Block S110B; back-coupling the second synchronization signal to generate a second self-receive signal in Block S120B; receiving the second self-receive signal in Block S130B; calculating a time-of-arrival of the second self-receive signal according to the second clock in Block S140B; receiving the first synchronization signal from the first node in Block S150B; and calculating a time-of-arrival of the first synchronization signal according to the second clock in Block S160B. The method S100 further includes calculating a time bias and a propagation delay between the pair of nodes based on the time-of-arrival of the first self-receive signal and the time-of-arrival of the second synchronization signal, the time-of-arrival of the second self-receive signal and the time-of-arrival of the first synchronization signal at the second node in Block S170.

As shown in FIG. 1B, a first variation of the method S100 includes: transmitting a coordination signal to a second node in the pair of nodes, the coordination signal indicating a second time to transmit a second synchronization signal in Block S102; transmitting a first synchronization signal at a first time according to a first clock of the first node in Block S110; back-coupling the first synchronization signal to generate a first self-receive signal in Block S120; receiving the first self-receive signal in Block S130; calculating a time-of-arrival of the first self-receive signal according to the first clock in Block S140; receiving the second synchronization signal, the second synchronization signal transmitted by the second node at the second time according to a second clock of the second node in Block S150; calculating a time-of-arrival of the second synchronization signal according to the first clock in Block S160; receiving a time-of-arrival of a second self-receive signal from the second node, the second self-receive signal back-coupled by the second node upon transmission of the second synchronization signal, the time-of-arrival of the second self-receive signal calculated according to the second clock in Block S142; receiving a time-of-arrival of the first synchronization signal from the second node, the time-of-arrival of the first synchronization signal calculated by the second node according to the second clock in Block S162; and calculating a first time bias between the first clock and the second clock and a first propagation delay between the first node and the second node based on the time-of-arrival of the first self-receive signal, the time-of-arrival of the first synchronization signal, the time-of-arrival of the second self-receive signal, and the time-of-arrival of the second synchronization signal in Block S170.

As shown in FIG. 1C, a second variation of the method S100 includes scheduling transmission of a first synchronization signal at a first time by a first node in Block S104; and scheduling transmission of a second synchronization signal at a second time by a second node in Block S106. The second variation of the method S100 also includes, after transmission of the first synchronization signal by the first node at the first time according to a first clock of the first node: receiving, from the first node, a time-of-arrival of a first self-receive signal according to the first clock, wherein the first self-receive signal is a back-coupled version of the first synchronization signal in Block S144; and receiving, from the second node, a time-of-arrival of the first synchronization according to a second clock of the second node in Block S164. The second variation of the method S100 further includes, after transmission of the second synchronization signal by the second node at the second time according to the second clock: receiving, from the second node, a time-of-arrival of a second self-receive signal according to the second clock, wherein the second self-receive signal is a back-coupled version of the second synchronization signal in Block S146; and receiving, from the first node, a time-of-arrival of the second synchronization signal according to the first clock in Block S166. The second variation of the method S100 also includes calculating a first time bias between the first clock and the second clock and a first propagation delay between the first node and the second node based on the time-of-arrival of the first self-receive signal, the time-of-arrival of the first synchronization signal, the time-of-arrival of the second self-receive signal, and the time-of-arrival of the second synchronization signal in Block S170.

A third variation of the method S100 includes: predicting a drift magnitude of the time bias between the first node and the second node between the first synchronization slot and the second synchronization slot based on a predictive drift model in Block S180; and setting a duration of the second synchronization slot equal to a sum of the maximum time bias uncertainty, the maximum accumulated drift magnitude over the duration of the first synchronization signal, the propagation delay, and a duration of the first synchronization signal in Block S190.

2. APPLICATIONS

Generally, the method S100 is executed by a pair of node devices (hereinafter “nodes”) in a network and/or remote servers in order to synchronize—such as to within one nanosecond—the clocks of the pair of nodes. The method S100 can calculate the time bias between clocks to within one nanosecond or with sub-nanosecond precision without prior information regarding the position of the two nodes (or even a requirement for the nodes to be stationary), without precise calibration of the hardware of either node, and with each node using standard electronic clock technology, such as a crystal oscillator clock. Additionally, the method S100 can utilize a frequency bandwidth of less than five megahertz. The small bandwidth of the method S100 enables transmission of frequency bands with high propagation range and/or penetration. Examples of such bands include frequency bands for cellular communications or public safety uses (e.g., the 902-928 MHz ISM band or sub-GHz cellular band). Thus, the method S100 enables applications such as time-synchronized distributed antenna systems, which can further enable precise localization of RF emitting devices, remote sensing of public safety infrastructure, and improvements in time-based data transfer protocols.

More specifically, the method S100 for synchronizing time between two nodes can be categorized as a two-way ranging and synchronization protocol. However, the method S100 is distinguished from other two-way ranging protocols in that each node transmits independently of other nodes at an absolute time (according to the clock of each node), as opposed to transmitting in response to receiving signals from a master node in the pair. Furthermore, each node is configured to execute a scheme to provide a local reference copy of the transmitted signal, within the transceiver hardware of the node, in order to provide a “time-of-departure” for each transmitted synchronization signal (via a back-coupling and/or reflection of a transmitted synchronization signal). The time-of-arrival of this reflected and/or back-coupled version of a transmitted synchronization signal (hereinafter a “self-receive signal”) includes the transmit and receive chain delays that would be incurred by a signal transmitted by another node, thereby enabling direct comparison of the time-of-arrival of various signals without precise characterization of these delays.

Blocks of the method S100 can be executed by a system, which can include a pair of node devices (i.e. nodes) in a mesh network, a master node and a slave node in a mesh network, and/or a remote server coordinating with a pair of nodes in a mesh network. Each node in the mesh network can include networking hardware, such as an antenna, transceiver hardware, an FPGA/DPS, a clock, and self-receive hardware (e.g. directional couplers, RF power splitters, combiners, circulators, etc.) further described below. However, the method S100 can also be performed between nodes in a wired network. The nodes in the mesh network are mutually connected to the internet or to a local area network such that the initial time bias between any pair of nodes in the mesh network is initially limited by the network time protocol (hereinafter “NTP”), or any other network time synchronization protocol. This time bias can range from tens of milliseconds to microseconds in many state-of-the-art networks.

Nodes in the mesh network communicate with each other using time division multiple access (hereinafter “TDMA”) and/or code division multiple access (hereinafter “CDMA”), thereby minimizing frequency bandwidth usage. Thus, the nodes can execute Blocks of the method S100 during one or more synchronization slots within a TDMA frame structure. In one implementation, synchronization slot duration and frame length are dynamically adjustable. In particular, the method S100 can leverage smaller initial time biases between a pair of nodes to decrease synchronization slot duration and further reduce uncertainty in the time bias calculation. Furthermore, the frame structure can include data transfer frames enabling the nodes to communicate with each other or with a remote server via the internet protocol suite in order to execute blocks of the method S100. For example, nodes can transfer a set of time-of-arrivals to a remote server for further processing according to the method S100.

Assuming an initial coarse clock synchronization between a pair of nodes (e.g. one to ten milliseconds), each of the pair of nodes transmits a synchronization signal to the other node at the beginning of the synchronization slot according to each node's clock. However, due to the relative time bias between each of the nodes, these synchronization signals are sent at times offset by the time bias between the two nodes. The synchronization signal is a frequency modulated, amplitude modulated, or phase modulated pseudorandom code or a combination of multiple codes (e.g., on multiple center carrier frequencies). Upon receiving the synchronization signal from the other node in the pair of nodes each node calculates a time-of-arrival (hereinafter “TOA”), such as by using the magnitude, time offset, and carrier phase of the autocorrelation peaks associated with one or more codes.

Upon sending the synchronization signal, each node includes self-receive hardware (shown in FIG. 3C) that back-couples and/or internally reflects a portion of the synchronization signal, thereby generating a self-receive signal, as shown in FIG. 3A. The self-receive hardware reflects and/or back-couples an attenuated repetition of the outgoing synchronization signal, which is then received by the same node that transmitted the synchronization signal. Each node then calculates a TOA for the self-receive signal after it has passed through the same receiver chain as an incoming signal would. Thus, the TOA of the self-receive signal acts as a time-of-departure of the synchronization signal offset by the receiver delay time of the transmitting node. The method S100 includes generating the self-receive signal such that the TOA of an incoming synchronization signal from another node can be directly compared to a time of departure of an outgoing node without the need for precise measurement and/or calibration of the receiver delay.

Each node in the pair then receives the synchronization signal from the other node and calculates a TOA for the synchronization signal. Therefore, by the end of the synchronization slot, each node has recorded a TOA for a self-receive signal and for a synchronization signal of the other node. These TOAs are then transmitted to one of the nodes or to a remote server to solve a system of two equations for two unknowns: the relative time bias between the pair of nodes (including receiver delay) and the propagation delay (i.e. transmission time) between the pair of nodes. This system of equations is solvable based on the reciprocity theorem of electromagnetism. The relative time bias between clocks of the two nodes can then be calculated and tracked by a master node or a remote server coordinating the mesh network. Alternatively, the relative time bias between nodes can then be reported to both nodes and one node (i.e. a slave node) can then synchronize its clock to match the other (i.e. a master node).

Although the first synchronization between the two node clocks according to the method S100 may be highly accurate, greater accuracy and reduced TDMA overhead is achieved through slot refinement, wherein the synchronization slot duration is reduced based on the smaller time bias between the two clocks. Therefore, in each successive synchronization slot for the pair of nodes the synchronization slot duration is decreased, thereby reducing sources of error such as accumulated jitter and/or environmental frequency drift that may occur during the synchronization process. For example, a typical quartz crystal oscillator may be expected to have a drift of six ppm which corresponds to six nanoseconds of drift during a millisecond duration synchronization slot. If the synchronization slot duration is instead reduced to 100 microseconds, the expected drift is reduced to less than one nanosecond.

By executing synchronization slot refinement techniques the system can also increase the duration of time within a frame for data transmission thereby improving data transfer rates for the nodes. Additionally, with a shorter synchronization slot duration, the method S100 can be performed at a higher rate of repetition thereby improving synchronization during changing environmental conditions.

After synchronization has occurred between a pair of nodes in the mesh network, the method S100 can be executed on a pairwise basis with other nodes in the mesh network in successive synchronization slots, thereby distributing time across the entire mesh network.

In some implementations, the method S100 leverages repetitive synchronization between two nodes and pairwise synchronization between multiple nodes in the network, in combination with environmental data recorded from inertial measurement units (hereinafter “IMUs”) or temperature sensors within each node, to create a predictive drift model for each node in the mesh network. The predictive drift model can characterize the time drift of a node as a function of measured environmental input data over time. The method S100 can also include triggering a synchronization slot in the next frame or altering the length of the synchronization slot based on the output of the predictive drift model. Alternatively, the method S100 can include adjusting a node's clock between synchronization slots.

3. EXAMPLES

The method S100 for synchronizing time between two node devices has wide ranging benefits including: more precise and flexible localization using a time synchronized mesh network of nodes; and improvements to time or frequency-based communication protocols.

In one example application of the method S100, a mesh network of nodes can perform the method S100 on a pairwise basis to synchronize their clocks across the mesh network. Each of the nodes in the mesh network can then individually calculate a TOA of a signal from an RF emitting source not included in the mesh network. By comparing TOAs between the nodes in the mesh network and performing multilateration, the location of the RF emitting source can be calculated to a degree of precision limited by temporal synchronization between nodes. For example, with temporal synchronization within one nanosecond, the location of an RF emitting source can be calculated to within 30 centimeters relative to the nodes in the mesh network. If one node in the mesh network acts as an anchor node and is calibrated with accurate global positional information, the absolute position of the RF emitting source can also be detected.

In one example application of mesh network-based localization, nodes are distributed throughout an area, within buildings, or even on cars driving on the streets of an urban area. The nodes execute Blocks of the method S100 to periodically synchronize their clocks and then collect TOA data for other RF emitting sources within the mesh area. The remote server coordinating with the nodes can subsequently transform these TOA data into three-dimensional geospatial locations of the RF emitting sources, thereby improving location services, asset tracking, and object detection. Thus, the remote server or master node can: calculate a relative location of nodes in the network based on the propagation delay between the nodes; receive, from several nodes in the mesh network, a time-of-arrival of a signal from a transmitting device; and calculate, via multilateration, a relative location estimate for the transmitting device accounting for the relative time bias between each of the nodes in the mesh network.

In another example, TDMA based protocols may be limited to slot durations of a few tens of milliseconds. Therefore, Blocks of the method S100 can be executed to synchronize clocks of devices within a wired or wireless network in order to more accurately calculate the time bias, enabling a reduction in dead time and interference between slots while improving bitrate and bandwidth due to reduced inefficiencies in the multiple access protocol.

4. NODE HARDWARE

As discussed above, the method S100 is executed by or via coordination with nodes in a network. Examples of hardware for wireless nodes are shown in FIGS. 3B and 3C. Generally, nodes include transmission components and receiver components, an FPGA or DSP configured to generate and process signals, a clock, and self-receive signal hardware.

In one implementation, the nodes transmit information wirelessly and as such include RF transceiver hardware such as a super heterodyne radio architecture and an Rx/Tx antenna shown in FIG. 3B. In this implementation, each node includes a “receive chain” and a “transmit chain.” The receive chain includes a pipeline of hardware components that process signals received from the Rx port. The transmit chain includes a pipeline of hardware components that process transmission signals generated by the FPGA or DSP and feed them into the Tx port. The receive chain and the transmit chain impart a “receive chain delay” and a “transmit chain delay” respectively. “Receive chain delay” and “transmit chain delay” describe the amount of time elapsed as a signal traverses either the receive chain or the transmit chain respectively.

In an alternate implementation, the nodes can communicate over a wired network. In this implementation, the nodes can include an I/O port and/or appropriate interface converters for communicating over any wired medium (e.g. ethernet, fiber optic, etc.) instead of an antenna. These interface converters can also impart a receive chain and transmit chain delays upon transmitted and received signals.

In one implementation, a node includes extant transceiver infrastructure such as cell tower, mobile phone, or any other RF transceiving device that has been adapted to execute the method S100. A cell tower or other extant transceiver can be adapted to execute the method S100 via software means. In alternate implementations, the nodes can include optimized hardware to improve the performance of the method S100, which can include impedance-matching networks at an antenna interface between a passive coupling device (e.g., coupling the transmit chain and the receive chain to the antenna).

In wireless node implementations, the FPGA or DSP of each node is configured to generate complex digital signals and output the generated signals to a DAC. The complex components of the digital signals represent the in-phase and quadrature portions (i.e. I/Q) of the analog signal to be generated by the DAC. Additionally, the FPGA or DSP of the node receives digital signals from the antenna of the node via an ADC and timestamps received synchronization signals according to an instant value of its own clock and a TOA calculation process further described below.

Each node includes a clock, such as a crystal oscillator clock or an atomic clock, which may be responsible for time-keeping functions at the node. The method S100 can be executed to synchronize the clocks of multiple nodes within a network. In particular, the method S100 can effectively synchronize crystal oscillator clocks that satisfy basic frequency stability, phase noise, and frequency requirements for wireless communication. In one example, the clock is a quartz crystal oscillator with an AT cut and a clock frequency of ten megahertz. However, a node can include a crystal oscillator of any frequency or cut assuming the aforementioned constraints are met.

Each node can include specific self-receive hardware, which back-couples and/or reflects a self-receive signal in Block S120. Generally, upon a node transmitting a synchronization signal, the self-receive hardware reflects and/or back-couples an attenuated repetition of the synchronization signal back to the receiving port of the transmitting node. By processing the transmitted synchronization signal via the receive chain and calculating a TOA for the self-receive signal, each node is able to timestamp transmitted signals delayed by the receive chain delay. When the node later receives a synchronization signal from another node, the TOA of the synchronization signal also incurs the same receive chain delay. Because the TOA of the self-receive signal and the TOA of the synchronization signal both include the receive chain delay they can be directly compared without precise calibration of the receive chain hardware.

In one implementation, the self-receive hardware includes an impedance mismatched directional coupler shown in FIG. 3C. The self-receive hardware can also include a variable impedance circuit, which may be software controlled to vary the gain of the reflected Tx signal into the Rx port. Depending on the specific hardware implementation of the node, similar impedance-matching can be applied to circulators, power splitters, or any other transmission line device, as further described below. In one implementation, a node can include standard antenna interface hardware with unintentionally mismatched antenna impedance (e.g., characterized by a voltage standing wave ratio greater than one), which can function as the self-receive hardware by reflecting synchronization signals at the antenna interface. Thus, each node need not include specialized self-receive hardware in order to execute the method S100.

The self-receive hardware defines a tuned reflection coefficient for signals incident to the interface between a passive coupling device and the antenna of the node, by including an impedance-matching network between the antenna and the passive coupling device. The impedance-matching network transforms the impedance of the antenna to an impedance that effects a precise reflection coefficient between the passive coupling device and the antenna. The reflection coefficient is selected such that the reflected power of the self-reflected signal is above the noise floor of the ADC and below the saturation voltage of the ADC.

In another implementation, the matching network can maintain the reflection coefficient for a wide range of signal frequencies. For example, the self-receive hardware can include a switchable wideband matching network to improve the consistency of the reflection coefficient across a wide range of transmit frequencies. In this example, the self-receive hardware can actively switch between multiple impedance-matching networks depending on the frequency of the signal being transmitted by the node.

In yet another implementation, the self-receive hardware can include a frequency multiplexer connected to multiple impedance-matching networks such that each distinct impedance-matching network encounters signals at a frequency that results in the precise reflection coefficient.

In an additional implementation, the system 100 can include an adaptive impedance-matching network that can adjust its impedance and therefore the reflection coefficient between the passive coupling device and the antenna. The node can then use the adaptive impedance-matching network to adjust the reflection coefficient ensuring reception of the self-reflected signal despite changes in the noise levels at the ADC due to changing gains in RF transmit chain and/or receive chain amplifiers and/or interfering signals in the frequency band or bands of interest.

In one implementation, the self-receive hardware includes a directional coupler as the passive coupling device. The directional coupler includes four ports, two for each of the coupled transmission lines. In this implementation, the transmit port and the antenna port are located on the same transmission line in the directional coupler while the receive port is located at the coupled port opposite the antenna, thereby receiving coupled power from the antenna and reflected power from the interface between the directional coupler and the antenna. Alternatively, the self-receive hardware can include a directional coupler wherein the receive port and the antenna port are located on the same transmission line. As a result, the receive port receives direct power from the antenna port and the reflected power from the antenna port. However, in this alternative implementation, the antenna port receives signals from the transmit port at a lower power due to coupling between the two transmission lines. Thus, a higher power at the transmit port results in the same power at the antenna port. However, in an implementation wherein the receive port and the antenna port are coupled to the same transmission line, the self-receive hardware can improve the sensitivity of the transceiver to signals received at the antenna.

In an alternative implementation, the self-receive hardware includes a power divider as the passive coupling device. The power divider divides power from an input port between two output ports. In this implementation, power input at the antenna port is split between the transmit port and the receive port, thus any power reflected at the interface between the antenna and the power divider is also divided between the transmit port and the receive port.

In yet another implementation, the self-receive hardware includes a circulator as the passive coupling device. The circulator can couple the transmit port directly to the antenna port while the antenna port is directly coupled to the receive port. Any reflections generated at the antenna port interface are then coupled back to the receive port.

In one implementation, each node can include a software-defined radio architecture performing the function of any of the hardware elements described above.

However, the self-receive hardware can include any software or hardware system for feeding a Tx signal into an Rx port of a node.

5. GROSS TIME SYNCHRONIZATION

Prior to execution of the method S100, the clocks of each node are coarsely synchronized using another time synchronization protocol. For example, the pair of nodes may be connected to the internet and in communication with one or more NTP servers. Thus, both node clocks may be synchronized to between one and ten milliseconds under some conditions. Alternatively, the nodes can communicate with a global navigation system (hereinafter “GNSS”) time synchronization server and may be synchronized to within one microsecond in some cases.

The method S100 leverages existing time synchronization protocols, such that the nodes are able to communicate with the same synchronization slot in a TDMA described below.

Thus, a remote server can: coarsely synchronize a first node's clock and a second node's clock; schedule transmission of a first synchronization signal within a first synchronization slot of a first frame, the first synchronization slot characterized by a first synchronization slot duration; schedule transmission of a second synchronization signal within a second synchronization slot of the first frame, the second synchronization slot characterized by the first synchronization slot duration; configure the second node to receive the first synchronization signal during the first synchronization slot of the first frame; and configure the first node to receive the second synchronization signal during the second synchronization slot of the first frame.

6. SLOT AND FRAME DEFINITION

Generally, in the method S100, communication between a pair of nodes occurs according to a slot and frame TDMA structure; examples of which are shown in FIGS. 4A and 4B. The TDMA structure can include one or more slots for both time synchronization via the method S100 and data transfer. In one implementation, shown in FIG. 4A, each TDMA frame includes a synchronization slot to synchronize each unique pair of nodes in the mesh network. For example, a mesh network including three nodes, n₁, n₂, and n₃, would include a first synchronization slot for the synchronization of n₁ and n₂; a second synchronization slot for n₁ and n₃; and a third synchronization slot for n₂ and n₃. Alternatively, as shown in FIG. 4B, the TDMA frame can include a synchronization slot for each node in the mesh network to transmit a synchronization signal to all other nodes in the mesh network. For example, in a mesh network including three nodes, n₁, n₂, and n₃, would include a first synchronization slot for n₁ to transmit a synchronization signal to n₂ and n₃; a second synchronization slot for n₂ to transmit a synchronization signal to n₁ and n₃; and a third synchronization slot for n₃ to transmit a synchronization signal to n₁ and n₂.

In one implementation, (e.g., if the initial synchronization accuracy achieved with NTP is sufficiently low), each synchronization slot is further divided into two subsequent sub-slots to reduce synchronization overhead. For example, in a synchronization for n₁ and n₂, n₁ transmits during the first sub-slot and n₂ receives the transmission from n₁. Then, in a second sub-slot, n₂ transmits and n₁ receives.

Each TDMA frame can also include a header indicating the presence and order of synchronization slots and data transfer slots for each node in the mesh network. In one implementation, each TDMA frame includes a set of synchronization slots. Alternatively, each TDMA frame may or may not include a set of synchronization slots according to the header for the TDMA frame. In yet another implementation, a TDMA frame can include multiple sets of synchronization slots. In one implementation, the header can include a coordination signal transmitted by a remote server or master node coordinating the TDMA protocol in order to communicate specific transmission times for synchronization signals to nodes in the mesh network.

The duration of synchronization slots and the buffer time between two synchronization slots corresponding to different pairs of nodes can vary depending on the implementation and factors such as the current known time bias between nodes and associated uncertainty. The method S100 can include adjusting synchronization slot duration based on the output of a predictive drift model, further described below. Furthermore, the synchronization slot duration can have a lower limit equal to a sum of the uncertainty in the time bias between the pair of nodes, the propagation time between the pair of nodes, the uncertainty of the propagation time between the pair of nodes, and the duration of the synchronization signal, of which the duration of the synchronization signal is typically the most significant factor. In one implementation the slot and the buffer between synchronization slots are both set to one millisecond. Additionally or alternatively, the method S100 can include adjusting the synchronization slot duration as a function of the signal-to-noise ratio of previously received synchronization signals in order to provide additional measurement acquisition time to determine TOA of signals with lower signal-to-noise ratios.

Additionally, the duration of each TDMA frame is subject to a set of practical constraints. The TDMA frame duration can have an upper limit defined by the expected drift of the least stable node clock in the node pair compared to the desired accuracy of the clock synchronization process. For example, if one node clock is known to drift one nanosecond per second due to accumulated jitter and the desired accuracy of the clock synchronization process is one nanosecond, the TDMA frame duration (or time between synchronization slots) would be limited at one second such that the time bias between nodes is not likely to exceed one nanosecond.

Furthermore, the TDMA frame duration has a lower limit based on the sum of the expected initial clock offset between the nodes, the total duration of the synchronization slots, the total duration of any data transfer slots, and the total duration of any time buffers between slots. The frame duration must be sufficiently long for each node to receive the initial synchronization signal from each of the other nodes in the network while accounting for a large initial time bias. Thus, if the nodes are initially synchronized using NTP, the frame duration should be at least tens of milliseconds to account for time biases that are typical between NTP synchronized clocks.

In one implementation, during initial synchronization of a set of nodes in a mesh network. The set of nodes are configured (e.g., by the remote server) to broadcast a coordination signal at a predetermined time according to the clock of each node. Given initial coarse time synchronization between nodes, the system can designate a master node as the first node to transmit a coordination signal to the other nodes. Each node in the mesh network can then confirm receipt of the master node's confirmation signal, thereby designating it as the master node. Thus, a first node in a pair of nodes can receive, from a second node in the pair of nodes, a confirmation of receipt of the coordination signal; and, in response to receiving confirmation of receipt of the coordination signal before the first node receives a second coordination signal from another node in the mesh network, designating the first node as a master node.

In an alternative implementation, the TOA of the coordination signal can be calculated at each node and sent to a remote server, which can then determine the master node in the mesh network. In yet another alternative implementation, the remote server designates the master node or coordinates the TDMA and the method S100 for the set of nodes in the mesh network.

When communicating within a TDMA synchronization slot, each node in a pair of nodes executing the method S100 can transmit synchronization signals and receive both synchronization signals from other nodes and self-receive signals during the synchronization slot. Therefore, in one implementation, the nodes execute a time division duplex (hereinafter “TDD”) or a half-duplex scheme to communicate within a synchronization slot. In this implementation, the synchronization signals transmitted by each node in a node pair are separated by a predetermined delay time (e.g., approximately half the synchronization slot duration) such that a first node in a pair of nodes can transmit a first synchronization signal and receive a self-receive signal corresponding to the first synchronization signal during a first division of the synchronization slot and subsequently receive a second synchronization signal (transmitted from a second node in the node pair) in a second division of the synchronization slot. In this implementation, nodes executing the method S100 can transmit a coordination signal including timing information such as slot division durations and timings to other nodes in the mesh network thereby establishing a half-duplex communication link between the two nodes.

In one implementation, the system can implement a full-duplex (i.e. frequency division duplex, hereinafter “FDD”) communication scheme whereby each node in a pair of nodes can transmit and receive synchronization signals simultaneously (e.g., wherein each synchronization signal is sent over multiple frequencies). Additionally, each node in a pair of nodes can be configured to simultaneously receive a self-receive signal back-coupled to the receiver of the node and a synchronization signal from a second node in the pair of nodes.

Therefore, in one implementation, a master node in a pair of nodes can transmit the coordination signal to the second node, wherein the coordination signal indicates a second time to transmit the second synchronization signal equal to a first time of transmission of the first synchronization signal; and establish a full-duplex communication link with the second node.

However, in applications wherein a larger bandwidth is available, the nodes can also communicate using frequency division multiple access or any other channel access method.

7. TIME BIAS AND PROPAGATION DELAY CHARACTERIZATION

As shown in FIG. 1A, each node in a pair of nodes can execute Blocks S110, S120, S130, S140, S150, and S160 to: transmit synchronization signals to each other at a predetermined transmission time that is evaluated according to the clock of each node; back-couple the synchronization signal via the self-receive hardware of the node to form a self-receive signal; and record TOAs for the back-coupled self-receive signal and the synchronization signal received from the opposite node in the pair of nodes. Thus, each node in a pair of nodes calculates two TOAs according to the node's clock. Each node can then send the TOA of the self-receive signal and the TOA of the synchronization signal (e.g., during a data transfer slot in the frame) either to one of the nodes in the pair of nodes or a remote server in order to calculate the time bias between the pair of nodes and the propagation delay between the pair of nodes, in Block S170.

Variations of the method S100 shown in FIGS. 1B and 1C include additional Blocks for coordinating a pair of nodes executing the method S100. In the first variation shown in FIG. 1B, a master node in the pair of nodes executes Blocks S102, S142, S162, and S170 in addition to the previously mentioned Blocks of the method S100 in order to: communicate a transmission time to the slave node; receive a TOA of the self-receive signal calculated by the slave node; receive a TOA of the synchronization signal calculated by the slave node; and calculate, at the master node, a time bias and a propagation delay between the master node and the slave node. Thus, the master node can communicate a transmission time to the slave node such that the master node can receive the synchronization signal transmitted by the slave node during a synchronization slot in the TDMA communication scheme.

In the second variation of the method S100 shown in FIG. 1C, a remote server or other computational device communicating with a pair of nodes executes Blocks of the method S104, S106, S144, S146, S164, S166, and S170 to: schedule a transmission time for each of the pair of nodes; receive TOAs from each of the nodes for a synchronization signal received by each of the pair of nodes, and a self-receive signal received at each of the pair of nodes; and calculate a time bias and a propagation delay between the pair of nodes.

Therefore, in any of the abovementioned variations, various entities included with the system can cooperate to execute a time synchronization protocol between a pair of nodes in a mesh network in order to calculate a time bias and a propagation delay between the pair of nodes. In one implementation, the time biases are tracked by a remote server and accounted for when managing a TDMA communication protocol, when executing multilateration for transmitting RF devices, or when executing any other protocol utilizing a time synchronized mesh network. Alternatively, one of the nodes can function as a master node and track time biases and propagation time between nodes in the mesh network.

In order to synchronize clocks of the pair of nodes relative to each other, each node in the node pair can concurrently or consecutively execute Blocks S110, S120, S130, S140, S150, and S160. For example, a first node in the pair of nodes can execute a first instance of these Blocks of the method S100 in S110A, S120A, S130A, S140A, S150A, and S160A, while a second node in the pair of nodes can execute a second instance of these Blocks in S110B, S120B, S130B, S140B, S150B, and S160B. For the purpose of clarity, description with reference to Blocks S110, S120, S130, S140, S150, or S160 can therefore apply to respective Blocks S110A, S120A, S130A, S140A, S150A, or S160A executed by the first node or respective Blocks S110B, S120B, S130B, S140B, S150B, or S160B executed by the second node.

Depending on the relative time bias between the two nodes, the first node can execute Blocks S110A, S120A, S130A, and S140A at any time relative to the second node executing Blocks S110B, S120B, S130B, and S140B (assuming both nodes are executing during the synchronization slot). However, the first node executes Blocks S150A and S160A to receive and calculate a TOA for a synchronization sent from the second node after the second node executes Blocks S110B, S120B, S130B, and S140B. Likewise, the second node executes Blocks S150B and S160B after the first node executes Blocks S110A, S120A, S130A, and S140A.

7.1 Transmission Times

In Blocks S110A and S110B each node in a pair of nodes transmits a synchronization signal at a predetermined or coordinated transmission time within a synchronization slot. In implementations wherein the nodes execute TDD, a first transmission time for Block S110A can be offset from a second transmission time for Block S110B by a transmission interval. Thus, each node in a pair of nodes transmits a synchronization signal at a predetermined transmission time according to each node's own clock. For example, if a first transmission time for a first node is set at 1:00:00, then the first node transmits the synchronization at 1:00:00 according to its own clock. Thus, from a third party perspective, each node transmits a synchronization signal offset by the time bias between the two nodes and offset by the transmission interval between the transmission times for each node.

Because each node transmits a synchronization signal at a predetermined transmission time, instead of transmitting a synchronization signal in response to receiving a signal from another node, such as is in typical two-way ranging protocols, the system executing the method S100 does not require precise characterization of receive chain delay, transmit chain delay, or intervening processing delay for each node in order to accurately calculate the time bias and propagation time between a pair of nodes.

7.2 Synchronization Protocol

Generally, each node generates a synchronization signal; concurrently or consecutively transmits the synchronization signal to the other node in a pair of nodes and receives a self-receive signal based on the transmitted synchronization signal; and receives a synchronization signal from the other node in the node pair. For ease of description, the nodes in a node pair will be referred to as n₁ and n₂. However, n₁ and n₂ can refer to any two-node pair in a mesh network (e.g., n_(i) and n_(j)). The method S100 includes calculating the relative time bias between n₁ and n₂ represented as (b₁−b₂) and the propagation delay τ. Based on the reciprocity theorem of electromagnetism, the propagation delay for a signal transmitted from n₁ to n₂ is equal to the propagation delay for a signal transmitted from n₂ to n₁ (τ_(1,2)=τ_(2,1)≡τ).

In calculating TOA such as in Blocks S140 and S160 the node can cross-correlate a received synchronization signal or self-receive signal with a template signal to determine a timestamp corresponding to the peak value of the autocorrelation function. Synchronization signals can include specifically chosen sequences that have a high autocorrelation value when the sequences are aligned and a low autocorrelation value otherwise. In one implementation, each node performs digital autocorrelation between received signals and template signals. Alternatively, each node performs analog autocorrelation between an analog conversion of the digital synchronization signal and the template signal. The latter accounts for the value of the samples in addition to the time durations between samples. Furthermore, the node can refine the calculated TOA based on phase information extracted from the synchronization signal. In one implementation, the nodes can transmit each synchronization signal over multiple frequency bands to increase the number of carrier phase measurements obtained by the receiving node and therefore increase the accuracy of the TOA calculation.

In Block S110A, at local time T₁ (i.e. a first transmission time), n₁ generates a digital synchronization signal at the FPGA or DSP. In one implementation, n₁ converts the digital synchronization signal to an analog synchronization signal for transmission via the antenna and propagates the converted synchronization signal along the transmit chain of n₁, thereby incurring a transmit chain delay, t₁. The analog signal then interacts with the self-receive hardware, which reflects or otherwise back-couples the synchronization signal into the Rx port of n₁ as shown in Block S120A. The self-receive signal (i.e. reflected or otherwise back-coupled) synchronization signal is substantially similar to the synchronization signal though the power of the two signals may differ. In one implementation, the power of the self-receive signal ensures the self-receive signal does not saturate the ADC of n₁.

In Block S130A, n₁ receives the self-receive signal as it is transmitting the synchronization signal. While the self-receive signal propagates through the receive chain of n₁ it incurs receive chain delay, r₁. In Block S140A, n₁ calculates a local TOA of its self-receive signal, S_(1,1), which is related to the unknown parameters of interest via the following equation:

S _(1,1) =T ₁ +t ₁ +r ₁.

Thus, S_(1,1) represents the TOA of the self-receive signal from n₁ according to the clock of n₁.

As shown in FIG. 5A, and in Blocks S110B, S120B, S130B, and S140B, n₂ executes the equivalent steps to those executed by n₁ in Blocks S110A, S120A, S130A, and S140A, thereby transmitting a second synchronization signal and calculating a TOA of the self-receive signal generated by n₂:

S _(2,2) =T ₂ +t ₂ +r ₂,

where T₂ is a local time at n₂ (i.e. the second transmission time), t₂ is that transmit chain delay of n₂, and r₂ is the receive chain delay of n₂. In implementations wherein the transmission times for each of the nodes is offset a transmission interval, Δ, T₂=T₁+Δ Thus, S_(2,2) represents the measured TOA of the self-receive signal from n₂ according to the clock of n₂.

As shown in FIG. 5B, and in Block S150A, n₁ receives the synchronization signal from n₂, which propagates through the receive chain of n₁ incurring receive chain delay, r₁. In Block S160A, n₁ then calculates a local TOA for the synchronization signal received from n₂, S_(1,2), which is expressed in the following equation:

S _(1,2) =T ₂ −b ₂ +t ₂ +τ+r ₁ +b ₁.

Thus, S_(1,2) represents the measured TOA of the synchronization signal from n₂ according to the clock of n₁.

In Blocks S150B and S160B, n₂ executes the equivalent steps to those executed by n₁ in Blocks S150A and S160A, thereby receiving a synchronization signal from n₁ and calculating a local TOA (at n₂) for the synchronization signal, S_(2,1), expressed as follows:

S _(2,1) =T ₁ −b ₁ +t ₁ +τ+r ₂ +b ₂.

Thus, S_(2,1) represents the measured TOA of the synchronization signal from n₁ according to the clock of n₂.

In Block S170, one of the nodes (e.g., a master node in the pair of nodes), either n₁ or n₂, or another computational device such as a remote server or a separate master node, collects the above TOAs, S_(1,1), S_(1,2), S_(2,2), and S_(2,1), to calculate τ and b₁−b₂+r₁−r₂, or the relative time bias plus the difference in receiver chain delay. Taking differences between the self-receive signal and the corresponding synchronization signal received at opposite nodes the following is derived:

S _(1,2) −S _(2,2) =τ+r ₁ +b ₁ −r ₂ −b ₂,

S _(2,1) −S _(1,1) =τ+r ₂ +b ₂ −r ₁ −b ₁.

Thus τ is calculated as:

τ=½[(S _(1,2) −S _(2,2))+(S _(2,1) −S _(1,1))]

and b₁−b₂+r₁−r₂ is calculated as:

b ₁ −b ₂ +r ₁ −r ₂=½[(S _(1,2) −S _(2,2))−(S _(2,1) −S _(1,1))].

Although a value for the pure relative time bias, b₁−b₂, cannot be calculated without the value of r₁−r₂, adjusting the clock of one node in the pair of nodes by b₁−b₂+r₁−r₂ ensures that any signal received by both nodes at the same time will receive the same timestamp at each node even if the instantaneous time at each clock differs by r₁−r₂ because the difference in receive chain delay will compensate for the post inaccuracy in the time bias. Alternatively, if r₁−r₂≈0, then:

b ₁ −b ₂≈½[(S _(1,2) −S _(2,2))−(S _(2,1) −S _(1,1))].

Once τ and b₁−b₂+r₁−r₂ have been calculated, the method S100 can also include synchronizing the time between n₁ and n₂ by adding b₁−b₂+r₁−r₂ to the clock of n₁ or n₂ to compensate for the initial time bias. In implementations with multiple nodes in the mesh network, one node is designated as a “master node” while the other nodes are designated “slave nodes.” Thus, in this implementation, the method S100 includes adjusting the clocks of the slave nodes to match the master node. Alternatively, a master node or remote server in the system can track the relative time bias of each node pair in a mesh network and compensate for the calculated time biases when executing processes that rely on precise time synchronization between nodes in the mesh network, such as time-based communication protocols or multilateration of other RF devices.

In one implementation, each node in the mesh network can average measurements to calculate a time bias and propagation delay when these variables are over-defined. For example, in a network including nodes n₁, n₂, and n₃, the time bias between n₁ and n₂ can also be calculated as follows:

b ₁ −b ₂≈½[(S _(1,3) −S _(3,3))−(S _(3,1) −S _(1,1))]−½[(S _(2,3) −S _(3,3))−(S _(3,2) −S _(2,2))].

Additionally, the method S100 can include indirectly calculating a time bias between nodes that do not have a direct communication line (e.g., due to an obstruction of communication between the two nodes). For example, in a network including nodes n₁, n₂, and n₃, the time bias between n₁ and n₃ can also be calculated as follows:

b ₁ −b ₃≈½[(S _(1,2) −S _(2,2))−(S _(2,1) −S _(1,1))]+[(S _(2,3) −S _(3,3))−(S _(3,2) −S _(2,2))].

Thus, n₂ and n₃ can execute another iteration of the method S100 and the system can calculate a time bias between the first clock of n₁ and the third clock of n₃ based on a sum of the first time bias (b₁−b₂) and the second time bias (b₂−b₃).

Furthermore, the method S100 can include calculating an instantaneous uncertainty in the relative time bias between nodes, δ(b₁−b₂) and an instantaneous uncertainty in the propagation delay δτ. Sources of uncertainty can include the propagated uncertainty based on the peak width of the autocorrelation function for each TOA calculation plus any expected phase noise that may occur at the clocks of each node between synchronization slots.

8. PHASE REFINEMENT

In addition to calculating the time bias and propagation time between nodes based on the synchronization protocol described above, the system can also refine the propagation time r by extracting carrier phase information from each of the four received signals, the first self-receive signal received at n₁, the first synchronization signal received at n₂ the second self-receive signal received at n₂, and the second synchronization signal received at n₁. Upon receiving the first self-receive signal n₁ detects a carrier phase of the first self-receive signal, which is represented as follows:

∠_(1,1)=(2πf _(c) T ₁+ϕ_(1,Rx))mod 2π,

where f_(c) is the carrier frequency and ϕ_(1,Rx) is the phase delay induced by the receive chain of n₁. Likewise, n₂ detects a carrier phase of the second self-receive signal, which is represented as follows:

∠_(2,2)=(2πf _(c) T ₂+ϕ_(2,Rx))mod 2π,

where f_(c) is the carrier frequency and ϕ_(1,Rx) is the phase delay induced by the receive chain of n₁. Additionally, n₁ can detect the carrier phase of the synchronization signal received from n₂, which is represented as follows:

∠_(1,2)=(2πf _(c)(T ₂+τ)+ϕ_(1,Rx))mod 2π.

Likewise, n₂ can detect the carrier phase of the synchronization signal received form n₁, which is represented as follows:

∠_(2,1)=(2πf _(c)(T ₁+τ)+ϕ_(2,Rx))mod 2π.

To refine the calculation of τ, the above relations can be rearranged as follows:

$\begin{matrix} {{\angle_{2,1} - \angle_{1,1}} = {\left( {{2\; \pi \; {f_{c}\left( {T_{1} + \tau} \right)}} + \varphi_{2,{Rx}} - {2\; \pi \; f_{c}T_{1}} - \varphi_{1,{Rx}}} \right)\; {mod}\; 2\; \pi}} \\ {{= {\left( {\varphi_{2,{Rx}} + {2\; \pi \; f_{c}\tau} - \varphi_{1,{Rx}}} \right)\; {mod}\; 2\; \pi}},} \end{matrix}$ $\begin{matrix} {{\angle_{1,2} - \angle_{2,2}} = {\left( {{2\; \pi \; {f_{c}\left( {T_{2} + \tau} \right)}} + \varphi_{1,{Rx}} + {2\; \pi \; f_{c}T_{2}} - \varphi_{2,{Rx}}} \right){mod}\; 2\; \pi}} \\ {= {\left( {\varphi_{1,{Rx}} + {2\; \pi \; f_{c\;}\tau} - \varphi_{2,{Rx}}} \right){mod}\; 2\; {\pi.}}} \end{matrix}$

Thus, the system can refine the value of τ, based on the synchronization protocol described above by calculating:

${{\left( {\angle_{2,1} - \angle_{1,1}} \right) + \left( {\angle_{1,2} - \angle_{2,2}} \right)} = {\left( {4\; \pi \; f_{c}\tau} \right)\; {mod}\; 2\; \pi}},{\tau \in {\left\{ {{\frac{\left( {\angle_{2,1} - \angle_{1,1}} \right) + \left( {\angle_{1,2} - \angle_{2,2}} \right) + {n\; 2\; \pi}}{4\; \pi \; f_{c}}\text{:}n} \in {\mathbb{Z}}} \right\}.}}$

Once several measurements are obtained for a range of carrier frequencies, f_(c), the system can statistically solve the set of equations (e.g., via a least squares calculation or other numerical methods) to refine the propagation time τ.

Thus, in one implementation, a remote server or master node executing Blocks of the method S100 can, after transmission of the synchronization signal by n₁ at T₁ according to a clock of n₁: receive, from n₁, a phase of a first self-receive signal according to n₁; and receive, from n₂, a phase of the synchronization signal from n₁ according to n₂; and, after transmission of the synchronization signal by n₂ at T₂ according to a clock of n₂: receive, from n₂, a phase of a second self-receive signal according to the clock of n₂; and receive, from n₁, a phase of the synchronization signal from n₂ according to the clock of n₁. The remote server or master node can then refine the first propagation delay based on the phase of the first self-receive signal, the phase of the first synchronization signal, the phase of the second self-receive signal, and the phase of the second synchronization signal.

9. CLOCK ADJUSTMENT

Once the system has executed the synchronization protocol and/or phase refinement for a pair of nodes, the system can shift the clock of one node in each pair of nodes by the calculated time bias in order to synchronize the clocks of the pair of nodes. Alternatively, as shown in FIG. 2 the system can maintain the time bias of each node relative to a master node or other time standard (e.g., the clock of a remote server) and correct scheduling times sent to each node and/or timestamps received from each node based on the latest calculated time bias of the node relative to the master node or time standard—such as by subtracting the time bias of the node from any timestamps calculated by the node. Furthermore, the system can predict a current time bias for each node based on the most recently calculated time bias and the predictive drift model for the node (further described below).

10. SYNCHRONIZATION SLOT REFINEMENT

In Block S180, the method S100 includes calculating a duration of a subsequent synchronization slot between previously synchronized nodes. Once a pair of nodes have executed the method S100, their clocks may be aligned to within one to ten nanoseconds. Thus, the synchronization slot can be made successively smaller as the synchronization between the two nodes improves.

The method S100 can also include calculating durations of subsequent synchronization slots based on a combination (e.g., a uncertainty propagated sum) of: the time bias uncertainty; the expected drift in the relative time bias between synchronization slots (calculated via the predictive drift model further described below); the propagation delay; the propagation delay uncertainty; the greater of the transmit chain delay of the first node or the transmit chain delay of the second node; the greater of the receive chain delay of the first node or the receive chain delay of the second node; the duration of the synchronization signal, and/or the expected change in the propagation delay due to relative movement between the pair of nodes.

In one implementation, the synchronization slot is additionally lengthened by a time buffer to ensure that each node receives a complete synchronization signal within the synchronization slot. When the synchronization slot duration is set equal to the above sum, the synchronization slot duration is long enough to allow a synchronization signal generated at one node to be fully received at the second node in a case where there is greater than expected drift or error in any of the aforementioned quantities.

In one example, receive chain delay and transmit chain delay may be negligible when compared to the synchronization signal duration. Therefore, the method S100 can include adding a time buffer of a duration approximating a maximum sum of a typical receive chain delay and a typical transmit chain delay to calculate a synchronization slot duration.

In one implementation, following a failure to synchronize caused by, for example, a greater than expected drift in the time bias between nodes or large increase in propagation delay, the method S100 can include inserting a synchronization slot in a subsequent TDMA frame and increasing the synchronization slot duration to provide more time to send and receive synchronization signals. The method S100 can increase the synchronization slot duration incrementally (e.g. in increments of ten microseconds) until the synchronization process is received or the method S100 can include extending the synchronization slot duration to a significantly longer duration to increase the likelihood of receiving a signal.

Thus, a remote server or master node executing Blocks of the method S100 can decrease the synchronization slot duration for a subsequent synchronization slot corresponding to a pair of nodes such that the updated synchronization slot duration is less than the initial synchronization slot duration but greater than a sum of a first uncertainty in the first time bias between the clocks of the pair of nodes and a synchronization signal duration. The remote server or master node can then schedule transmission times for each node in the pair of nodes within the updated synchronization slot duration.

Furthermore, the system can schedule the transmission times for synchronization signals in a pair of nodes based on a recently calculated time bias between the pair of nodes. Thus, the system can account for the relative time bias between nodes in order to ensure that both nodes are transmitting during the next synchronization slot despite the reduced duration of the synchronization slot.

11. PREDICTIVE DRIFT MODEL

In one variation of the method S100, the master node or remote server executing Block S170, can also execute a predictive drift model based on successive time bias and propagation delay characterization between nodes in the mesh network. Generally, the predictive drift model characterizes the drift of a node's clock (e.g. a crystal oscillator clock) relative to the clocks of other nodes in the network or relative to a time standard (e.g., UTC), as a function of environmental factors such as temperature, humidity, movement, and vibration. More specifically, the method S100 can include recording environmental data at each node and observing relative drift calculated over multiple prior synchronization slots to train a predictive drift model of the drift of particular nodes in the mesh network. The predictive drift model outputs the expected drift of one node in the mesh network relative to a master node or a time at a remote server based on the time elapsed since the most recent synchronization slot for the node and recent environmental data collected at the node.

The system can maintain predictive drift models characterizing the drift of each pair of nodes in a mesh network relative to each other. Alternatively, the master node and/or remote server can compress the pairwise predictive drift models and instead maintain a single predictive drift model for each node in the mesh network relative to the clock at a master node in the mesh network or a time standard accessed at a remote server.

In one implementation, the method S100 includes periodically adjusting a node's clock time between synchronization cycles according to the output of the predictive drift model (e.g. updating the node's clock value every millisecond according to the predicted drift at the node). Alternatively, the system can calculate an estimate for the relative time bias of a node between synchronization cycles based on the predictive drift model of the node.

Furthermore, the method S100 can include adjusting the duration of a subsequent synchronization slot according to the output of the predictive drift model. Furthermore, the method S100 can include triggering the inclusion of a synchronization slot in a TDMA frame in response to a predicted drift value being greater than a predefined drift threshold (e.g. triggering the inclusion of a synchronization slot in the TDMA frame when the predicted drift model predicts greater than one nanosecond of drift between the slave node clock and the master node clock).

The predictive drift model can be a combination of a set of physical models representing pertinent environmental parameters. In one implementation, the predictive drift model includes a temperature model for drift as a function of temperature and crystal cut. The temperature model can include a set of approximate models, each model relating drift (e.g. in parts-per-million) to temperature for a corresponding crystal cut. In one implementation, the method S100 can include classifying the crystal oscillator of each node according to observed temperature drift over multiple synchronization slots to select a polynomial approximation of temperature drift associated with a particular node as a function of time.

The predictive drift model can also include predictive models for thermal hysteresis, ambient pressure, humidity, electric and/or magnetic field strength, drive level for the crystal oscillator, and/or reference voltage for the crystal oscillator.

In one implementation, the predictive drift model can output a confidence interval for expected drift based on randomly distributed phase noise in a crystal oscillator and systematic drift. The predictive drift model can calculate a distribution of phase noise for the crystal oscillator based on factors such as reference source noise, power supply noise, vibration induced noise, and/or acceleration induced noise.

Thus, the system can execute successive iterations of the method S100 for a node in the mesh network and record a time series of relative time biases between the node and a reference time, which can be a time maintained at a master node or a remote server. Concurrently, the system can record a time series of environmental data from the node. The system can then correlate the time series of environmental data with the time series of time biases to predict the drift of the time bias of the node based on changes in the environmental data. In this manner, the system can calculate a temperature correlation between the temperature at a particular node and the drift rate exhibited by the node based on a time series of temperature data and a time series of time biases of the node. Additionally or alternatively, the system can calculate a movement correlation between the movement of a node (e.g., acceleration measured by an IMU at the node) and the drift rate exhibited by the node based on a time series of movement data and a time series of time biases of the node.

In yet another implementation, the system can detect the movement of a node by recording a time series of the propagation delays between one node in a mesh network and other nodes of the mesh network. Thus, the system can: calculate a first relative location of a node based on a first propagation delay calculated during a first iteration of the method S100; calculate a second relative location of a node based on a second propagation delay calculated during a second iteration of the method S100; and calculate a movement correlation of drift of the clock of the node based on the first relative location, the second relative location, the first time bias, and the second time bias

In one implementation, the system can: calculate the time bias uncertainty for a recently calculated time bias of a first node relative to a second node; predict the drift magnitude of the time bias at the node based on a predictive drift model; predict a change in the propagation delay between the first node and the second node based on motion data from the first node and the second node; and set a synchronization slot duration equal to a sum of the time bias uncertainty plus the drift magnitude plus the propagation delay plus the change in the propagation delay. Therefore, the system can ensure reception of subsequent synchronization signals transmitted between a node pair by accounting for relative drift between a pair of nodes in the mesh network.

12. RESYNCHRONIZATION TRIGGERS

One variation of the method S100 includes, triggering the inclusion of one or more synchronization slots within a subsequent TDMA frame based on the output of the predictive drift model or the detection of an environmental change. In one implementation, the method S100 includes triggering the inclusion of a synchronization slot upon the predictive drift model outputting a predicted drift magnitude greater than a threshold drift value (e.g. greater than ten nanoseconds of drift since the last synchronization slot). Alternatively, the method S100 can include triggering the inclusion of a synchronization slot based on observed temperature change (e.g. via a digital thermometer at the node) and/or acceleration or vibration data (e.g. via IMU data collected at the node). Furthermore, the method S100 can include triggering the inclusion of a synchronization slot based on movement of one node relative to another. The method S100 can include detecting movement via an IMU at the node, a measurement of the doppler shift for incoming signals from the node, or via multilateration between the nodes.

In one example, the system can: after an initial iteration of the method S100, measure a temperature at a first node in the node pair; calculate a time bias drift based on the temperature and a temperature drift model; and schedule a synchronization slot in the next TDMA frame based on the time bias drift.

Generally, when method S100 triggers the inclusion of a synchronization slot, the synchronization slot is included in the subsequent TDMA frame and indicated in the header of the TDMA frame.

The systems and methods described herein can be embodied and/or implemented at least in part as a machine configured to receive a computer-readable medium storing computer-readable instructions. The instructions can be executed by computer-executable components integrated with the application, applet, host, server, network, website, communication service, communication interface, hardware/firmware/software elements of a user computer or mobile device, wristband, smartphone, or any suitable combination thereof. Other systems and methods of the embodiment can be embodied and/or implemented at least in part as a machine configured to receive a computer-readable medium storing computer-readable instructions. The instructions can be executed by computer-executable components integrated by computer-executable components integrated with apparatuses and networks of the type described above. The computer-readable medium can be stored on any suitable computer readable media such as RAMs, ROMs, flash memory, EEPROMs, optical devices (CD or DVD), hard drives, floppy drives, or any suitable device. The computer-executable component can be a processor but any suitable dedicated hardware device can (alternatively or additionally) execute the instructions.

As a person skilled in the art will recognize from the previous detailed description and from the figures and claims, modifications and changes can be made to the embodiments of the invention without departing from the scope of this invention as defined in the following claims. 

I claim:
 1. A method for characterizing time bias and propagation delay comprising: scheduling transmission of a first synchronization signal at a first time by a first node; scheduling transmission of a second synchronization signal at a second time by a second node; after transmission of the first synchronization signal by the first node at the first time according to a first clock of the first node: receiving, from the first node, a time-of-arrival of a first self-receive signal according to the first clock, wherein the first self-receive signal is a back-coupled version of the first synchronization signal; and receiving, from the second node, a time-of-arrival of the first synchronization according to a second clock of the second node; after transmission of the second synchronization signal by the second node at the second time according to the second clock: receiving, from the second node, a time-of-arrival of a second self-receive signal according to the second clock, wherein the second self-receive signal is a back-coupled version of the second synchronization signal; and receiving, from the first node, a time-of-arrival of the second synchronization signal according to the first clock; and calculating a first time bias between the first clock and the second clock and a first propagation delay between the first node and the second node based on the time-of-arrival of the first self-receive signal, the time-of-arrival of the first synchronization signal, the time-of-arrival of the second self-receive signal, and the time-of-arrival of the second synchronization signal.
 2. The method of claim 1, further comprising: calculating a relative location of the first node and a relative location of the second node based on the first propagation delay; receiving, from the first node, a first time-of-arrival of a signal from a transmitting device; receiving, from the second node, a second time-of-arrival of the signal from the transmitting device; and estimating a relative location of the transmitting device based on a difference between the first time-of-arrival of the signal and the second time-of-arrival of the signal.
 3. The method of claim 1, further comprising: calculating a first relative location of the second node based on the first propagation delay; scheduling transmission of a third synchronization signal by the first node; scheduling transmission of a fourth synchronization signal by the second node; after transmission of the third synchronization signal by the first node at the third time according to the first clock: receiving, from the first node, a time-of-arrival of a third self-receive signal according to the first clock, wherein the third self-receive signal is a back-coupled version of the third synchronization signal; and receiving, from the second node, a time-of-arrival of the third synchronization according to a second clock of the second node; after transmission of the fourth synchronization signal by the second node at the fourth time according to the second clock: receiving, from the second node, a time-of-arrival of a fourth self-receive signal according to the second clock, wherein the fourth self-receive signal is a back-coupled version of the fourth synchronization signal; and receiving, from the first node, a time-of-arrival of the fourth synchronization signal according to the first clock; based on the time-of-arrival of the third self-receive signal, the time-of-arrival of the third synchronization signal, the time-of-arrival of the fourth self-receive signal, and the time-of-arrival of the fourth synchronization signal, calculating a second time bias between the first clock and the second clock and a second propagation delay between the first node and the second node; and calculating a second relative location of the second node based on the second propagation delay; and calculating a movement correlation of drift of the second clock relative to the first clock based on the first relative location, the second relative location, the first time bias, and the second time bias.
 4. The method of claim 1: further comprising coarsely synchronizing the first clock and the second clock; wherein scheduling transmission of the first synchronization signal further comprises scheduling transmission of the first synchronization signal within a first slot of a first frame, the first slot characterized by a first slot duration; wherein scheduling transmission of the second synchronization signal further comprises scheduling transmission of the second synchronization signal within a second slot of the first frame, the second slot characterized by the first slot duration; further comprising configuring the second node to receive the first synchronization signal during the first slot of the first frame; and further comprising configuring the first node to receive the second synchronization signal during the second slot of the first frame.
 5. The method of claim 4, further comprising: scheduling transmission of a third synchronization signal by the first node at a third time within a first slot of a second frame, the third time accounting for the first time bias the first slot of the second frame characterized by a second slot duration, the second slot duration: greater than a sum of a first uncertainty in the first time bias between the first clock and the second clock and a synchronization signal duration; and less than the first slot duration; scheduling transmission of a fourth synchronization signal by the second node at a fourth time within a second slot of the second frame, the fourth time accounting for the first time bias the second slot of the second frame characterized by the second slot duration; configuring the second node to receive the third synchronization signal during the first slot of the second frame; configuring the first node to receive the fourth synchronization signal during the second slot of the second frame; and after transmission of the third synchronization signal by the first node at the third time according to the first clock: receiving, from the first node, a time-of-arrival of a third self-receive signal according to the first clock, wherein the third self-receive signal is a back-coupled version of the third synchronization signal; and receiving, from the second node, a time-of-arrival of the third synchronization according to the second clock; after transmission of the fourth synchronization signal by the second node at the fourth time according to the second clock: receiving, from the second node, a time-of-arrival of a fourth self-receive signal according to the second clock, wherein the fourth self-receive signal is a back-coupled version of the fourth synchronization signal; and receiving, from the first node, a time-of-arrival of the fourth synchronization according to the first clock; and based on the time-of-arrival of the third self-receive signal, the time-of-arrival of the third synchronization signal, the time-of-arrival of the fourth self-receive signal, and the time-of-arrival of the fourth synchronization signal, calculating a second time bias between the first clock and the second clock and a second propagation delay between the first node and the second node, the second time bias characterized by a second uncertainty less than the first uncertainty.
 6. The method of claim 1, further comprises calculating the first time bias between the first clock and the second clock and the first propagation delay between the first node and the second node based on the first time and the second time according to an ideal clock.
 7. The method of claim 1, further comprising: scheduling transmission of a third synchronization signal by the first node; scheduling transmission of a fourth synchronization signal by the second node; after transmission of the third synchronization signal by the first node at the third time according to the first clock: receiving, from the first node, a time-of-arrival of a third self-receive signal according to the first clock, wherein the third self-receive signal is a back-coupled version of the third synchronization signal; and receiving, from the second node, a time-of-arrival of the third synchronization according to a second clock of the second node; after transmission of the fourth synchronization signal by the second node at the fourth time according to the second clock: receiving, from the second node, a time-of-arrival of a fourth self-receive signal according to the second clock, wherein the fourth self-receive signal is a back-coupled version of the fourth synchronization signal; and receiving, from the first node, a time-of-arrival of the fourth synchronization signal according to the first clock; based on the time-of-arrival of the third self-receive signal, the time-of-arrival of the third synchronization signal, the time-of-arrival of the fourth self-receive signal, and the time-of-arrival of the fourth synchronization signal, calculating a second time bias between the first clock and the second clock and a second propagation delay between the first node and the second node; and based on the first time bias and the second time bias characterizing a first drift of the second clock relative to the first clock.
 8. The method of claim 7, further comprising: scheduling transmission of a fifth synchronization signal by the first node and scheduling transmission of a sixth synchronization signal by the second node based on the first drift.
 9. The method of claim 7, further comprising: receiving from the second node a first temperature of the second clock, the second temperature recorded at the second time; receiving from the second node a second temperature of the second clock, the fourth temperature recorded at the fourth time; and calculating a temperature correlation of drift of the second clock relative to the first clock based on the first temperature, the second temperature, and the first drift.
 10. The method of claim 9, further comprising: receiving from the second node a third temperature of the second clock, the third temperature recorded after the fourth time; calculating a second drift based on the third temperature and the temperature correlation of drift of the second clock relative to the first clock; scheduling transmission of a fifth synchronization signal by the first node and scheduling transmission of a sixth synchronization signal by the second node based on the second drift.
 11. The method of claim 1, further comprising: after transmission of the first synchronization signal by the first node at the first time according to a first clock of the first node: receiving, from the first node, a phase of a first self-receive signal according to the first node; and receiving, from the second node, a phase of the first synchronization signal according to the second node; after transmission of the second synchronization signal by the second node at the second time according to the second clock: receiving, from the second node, a phase of a second self-receive signal according to the second clock; and receiving, from the first node, a phase of the second synchronization signal according to the first clock; and refining the first time bias and the first propagation delay based on the phase of the first self-receive signal, the phase of the first synchronization signal, the phase of the second self-receive signal, and the phase of the second synchronization signal.
 12. The method of claim 1, further comprising: scheduling transmission of a third synchronization signal at a third time by a third node; after transmission of the second synchronization signal by the second node at the second time according to the second clock, receiving, from the third node, a second time-of-arrival of the second synchronization signal according to a third clock of the third node; after transmission of the third synchronization signal by the third node at the third time according to the third clock: receiving, from the third node, a time-of-arrival of a third self-receive signal according to the third clock, wherein the third self-receive signal is a back-coupled version of the third synchronization signal; and receiving, from the second node, a time-of-arrival of the third synchronization signal according to the second clock; calculating a second time bias between the second clock and the third clock and a second propagation delay between the second node and the third node based on the time-of-arrival of the second self-receive signal, the second time-of-arrival of the second synchronization signal according to the third clock, the time-of-arrival of the third self-receive signal, and the time-of-arrival of the third synchronization signal; and calculating a third time bias between the first clock and the third clock based on a sum of the first time bias and the second time bias.
 13. A method for calculating time bias between a pair of nodes in a network comprising, at a first node in the pair of nodes: transmitting a coordination signal to a second node in the pair of nodes, the coordination signal indicating a second time to transmit a second synchronization signal; transmitting a first synchronization signal at a first time according to a first clock of the first node; back-coupling the first synchronization signal to generate a first self-receive signal; receiving the first self-receive signal; calculating a time-of-arrival of the first self-receive signal according to the first clock; receiving the second synchronization signal, the second synchronization signal transmitted by the second node according at the second time according to a second clock of the second node; calculating a time-of-arrival of the second synchronization signal according to the first clock; receiving a time-of-arrival of the first synchronization signal from the second node, the time-of-arrival of the first synchronization signal calculated by the second node according to the second clock; receiving a time-of-arrival of a second self-receive signal from the second node, the second self-receive signal back-coupled by the second node upon transmission of the second synchronization signal, the time-of-arrival of the second self-receive signal calculated according to the second clock; and calculating a first time bias between the first clock and the second clock and a first propagation delay between the first node and the second node based on the time-of-arrival of the first self-receive signal, the time-of-arrival of the first synchronization signal, the time-of-arrival of the second self-receive signal, and the time-of-arrival of the second synchronization signal.
 14. The method of claim 13: wherein transmitting the coordination signal to the second node further comprises transmitting the coordination signal to the second node, the coordination signal indicating the second time to transmit the second synchronization signal, the second time offset from the first time by an transmission interval when evaluated by an ideal clock; and establishing a half-duplex communication link with the second node.
 15. The method of claim 13, further comprising: receiving, from the second node confirmation of receipt of the coordination signal; and in response to receiving the confirmation of receipt of the coordination signal before receiving a second coordination signal from the second node, designating the first node as a master node.
 16. The method of claim 13: wherein transmitting the coordination signal to the second node further comprises transmitting the coordination signal to the second node, the coordination signal indicating the second time to transmit the second synchronization signal, the second time equal to the first time when evaluated by an ideal clock; and further comprising establishing a full-duplex communication link with the second node.
 17. A method for calculating time bias between a pair of nodes in a network, including: at a first node in the pair of nodes, during a first synchronization slot: transmitting a first synchronization signal at a first time according to a first clock of the first node; back-coupling the first synchronization signal to generate a first self-receive signal; receiving the first self-receive signal; calculating a time-of-arrival of the first self-receive signal according to the first clock; receiving a second synchronization signal from the second node; and calculating a time-of-arrival of the second synchronization signal according to the first clock; at the second node, during the first synchronization slot: transmitting the second synchronization signal at a second time according to a second clock of the second node; back-coupling the second synchronization signal to generate a second self-receive signal; receiving the second self-receive signal; calculating a time-of-arrival of the second self-receive signal according to the second clock; receiving the first synchronization signal from the first node; and calculating a time-of-arrival of the first synchronization signal according to the second clock; calculating a time bias and a propagation delay between the pair of nodes based on the time-of-arrival of the first self-receive signal and the time-of-arrival of the second synchronization signal, the time-of-arrival of the second self-receive signal and the time-of-arrival of the first synchronization signal at the second node.
 18. The method of claim 17, wherein calculating the time bias and the propagation delay between the pair of nodes further comprises: calculating the time bias by solving for b₁−b₂ in: b ₁ −b ₂=½[(S _(1,2) −S _(2,2))−(S _(2,1) −S _(1,1))] calculating the propagation delay by solving for r in: τ=½[(S _(1,2) −S _(2,2))+(S _(2,1) −S _(1,1))] wherein S_(1,1) represents the time-of-arrival of the first self-receive signal, S_(1,2) represents the time-of-arrival of the second synchronization signal, S_(2,1) represents the time-of-arrival of the first synchronization signal, and S_(2,2) represents the time-of-arrival of the second self-receive signal.
 19. The method of claim 17, further comprising: calculating a time bias uncertainty between the first node and the second node; predicting a drift magnitude of the time bias between the first node and the second node between the first synchronization slot and a second synchronization slot based on a predictive drift model; predicting a change in the propagation delay between the first node and the second node based on a motion data from the first node and the second node; and setting a synchronization slot duration equal to a sum of the time bias uncertainty plus the drift magnitude plus the propagation delay plus the change in the propagation delay.
 20. The method of claim 19, further comprising: at the first node, during a second synchronization slot characterized by the synchronization slot duration: transmitting a third synchronization signal at a third time according to the first clock; back-coupling the third synchronization signal to generate a third self-receive signal; receiving the third self-receive signal; calculating a time-of-arrival of the third self-receive signal according to the first clock; receiving a fourth synchronization signal from the second node; and calculating a time-of-arrival of the fourth synchronization signal according to the first clock; at the second node, during the second synchronization slot: transmitting the fourth synchronization signal at a fourth time according to the second clock; back-coupling the fourth synchronization signal to generate a fourth self-receive signal; receiving the fourth self-receive signal; calculating a time-of-arrival of the fourth self-receive signal according to the second clock; receiving the third synchronization signal from the first node; and calculating a time-of-arrival of the third synchronization signal according to the second clock; calculating a second time bias and a second propagation delay between the pair of nodes based on the time-of-arrival of the third self-receive signal, the time-of-arrival of the fourth synchronization signal, the time-of-arrival of the fourth self-receive signal, and the time-of-arrival of the third synchronization signal at the second node. 